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ABSTRACT 

The measurement of channel gains is a widely used method for de- 
termining properties of substances. A common way of performing 
this measurement is by emitting a harmonic signa] with known am- 
plitude followed by a filtering of the received signal, which yields 
an estimation of the channel gain, and hence a quantification of 
an unknown substance property. These types of signals are unfor- 
tunately sensitive to various kinds of noise, making robust imple- 
mentations difficult. The wavelet transform, being fast and highly 
adaptable, is proposed as a tool for making these estimation ro- 
bust. This paper demonstrates that by combining a series of care- 
fully designed signals and the versatility of the wavelet transform, 
it is possible, with low computational complexity, to make robust 
estimations of channel gains. The concern of this paper is mea- 
surement of channel gains, but the method might very well apply 
to other areas. 



1. INTRODUCTION 

One of the ways of determining the density, transparency or thick- 
ness of a substance or material, such as smoke, wine, glass, paper, 
plastic is channel gain measurements; by transmitting a signal with 
a known intensity through the substance or material, and estimat- 
ing the resulting intensity, it is possible to determine the density, 
transparency or thickness. Channel gain measurements are also 
used for ascertaining the surface properties of, distance to, or mere 
presence of objects. In this case an emitted signal is reflected onto 
a receiver. An example is a system for automatically opening a 
door whenever a person is present in from of it. 

1.1. The previous solutions 

A typical way of making this type of measurements is emitting 
a simple signal, such as a harmonic signal, since this is easily 
constructed with analog electronics. Moreover the intensity of a 
harmonic signal is found by a plain band pass filtering. Alterna- 
tively a non-structured signal, like radiation or light with constant 
intensity, can be emitted. The estimation of the received inten- 
sity is in this case particularly simple. However, the use of simple 
signals makes the channel gain measurement sensitive to the envi- 
ronment, since most electrical apparatuses found in industry and 
at home emits (often unintentionally) all kinds of simple signals, 
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both visually and audibly. For instance a remote control produces 
as series of infra red signal, while a television set emits a rather in- 
tense, high frequent sound. For every new application that utilizes 
channel gain measurements in this fashion it is therefore necessary 
with a great deal of testing and fine-tuning. A necessity which is 
both costly and time consuming. 



1.2. The proposed solution 

Since the simplicity of the signals is the Achilles' heel with re- 
spect to robustness of the estimations of the channel gain, a logical 
approach is to introduce more complex signals. That is signals 
carefully designed such thai they are easily recognized even when 
overlaid with severe noise. Adaptive design of complex signals is 
most easily done digitally, and the proposed solution therefore as- 
sumes the possibility of digital signal processing. The key question 
then becomes what signal is best suited for transmission under the 
given conditions, including the overlaid noise, the characteristics 
of the electrical components, the properties of the transparent or 
reflecting material, and the capacity of the signal processor. The 
noise is a priori unknown, although in many applications some 
types of noise are likely to occur. Due to the low cost of the elec- 
trical components the characteristics of these can vary. quite a lot, 
even for two seemingly identical components. The low cost pri- 
ority also result in a fairly low signal processor capacity. In some 
applications the transparent or reflecting material is well-known, 
in others it is unknown. The large uncertainty on some of the main 
factors makes it virtually impossible to design 'the best' signal. 
For not only should the signal be immune to many types of noise, 
it is also subject to a trade off between recognizability (which often 
implies complexity), robustness, and processor capacity. 

Instead of transmitting a simple signal, a signal designed to 
have certain properties are 'protected' by an inverse transform be- 
fore transmission. The received signal is transformed resulting in 
the original signal overlaid with noise. Since the wavelet transform 
can produce predefined trade offs between time and frequency in- 
formation, and since the design of the original signal is completely 
free, this method can be adapted to virtually any type of noise. 
Because the transform is linear, is energy preserving, and has per- 
fect reconstruction, and because the original signal is known, it is 
easy to determine a number of properties of the current overlaid 
noise. This information can be used for automated, online adap- 
tation. Good introductions to the wavelet theory are (10] and [9j. 
A more rigorous, mathematical treatment of the subject is given in 



2. THE IDEA 

A digital signal s 0 is generated, processed and transmitted from 
emitter to receiver. The results is another signal s r , which is also 
processed to determine the intensity of transmitted signal, and hence 
the channel gain. The wavelet transform included in the processing 
is the discrete wavelet packet transform. The result of transform- 
ing, denoted by W(), is the representationof the signals in some 
basis, and the inverse transform, denoted W( ), is the reconstruc- 
tion of the signal from whatever basis it is represented in. 

The starting point is the signal so> the structure of which even- 
tually will determine the quality of the measurement, and from this 
the signal to be transmitted s t is constructed by an inverse wavelet 
transform (from whatever basis the signal so is represented in) fol- 
lowed by an affine mapping to adjust it to the range of the emitter 

s t = aVV(so) +/91. 

The transmission signal is given by s r = X(s t ), where T is the 
transfer function from emitter to receiver, including the character- 
istics of both components. T is assumed to be a constant transfer 
function, but the approach can easily be extended to handling dy- 
namical ones as well. Assume that the transmission dampens the 
signal and adds noise, that is T(x) = Gx + e t . Then the wavelet 
transform of the received signal becomes 

W(s r ) = W(C(aVV(s 0 ) + 01) + e t ) 
= G(as 0 +/?W(1)) + W(eO, 

where G is the channel gain. Since a and 0 are known this reduces 
to 



— Gs n + e w 



(1) 



where both s r (hence also s wr ) and s„ are known signals. This 
equation has three major degrees of freedom; the choice of original 
signal, the choice of wavelet transform, and the choice of solution 
method. The goal is a good estimation of the channel gain G, and 
by exploiting all three degrees of freedom this is possible even for 
severe noise conditions. 

2.1. The wavelet transform 

The wavelet transform has two purposes in this framework: Pro- 
ducing the emitted signal and transforming the noise contribution. 
The former is determinant for the energy consumption of the emit- 
ter, and for the emitters influence on the surroundings, while the 
latter has direct influence on the quality of the estimation of G. 
In this paper we will be concerned solely with the quality of the 
estimation. 

The wavelet transform itself has a number of degrees of free- 
dom, which can be exploited in the adaptation. Apart from the 
obvious choice of filter and choice of basis, the transform type is 
important; an integer-to-integer transform [1] might be attractive 
since the signals are digital (this would eliminate the quantization 
errors), the lifting technique [3, 8] could prove superior in adapt- 
ing the transform to signal and/or noise, while a transform on finite 
fields [5] perhaps is the best solution to the dynamic range problem 
occurring in fixed point arithmetic. Another transform property is 
the handling of the ends of the signal. This is important since W 
is a wavelet packet transform, which by repeated transformations 
yields several elements, each with two ends. The 'ends* problem is 



traditionally dealt with either by periodization or by mirroring. Al- 
ternatively boundary filters [6] or time-varying transform [7] could 
be used to reduce the influence of the noise at the ends. 

The choice of filter is usually be restricted to orthogonal filters, 
since non-orthogonal filters tend to have asymmetric frequency 
responses rendering the result of repeated transforms inhomoge- 
neous with respect to gain in the various frequency bands. This 
property might be useful for suppressing noise, however, but this 
has, to the best of the authors 1 knowledge, not yet been investi- 
gated. 

2.2. Solving the equation 

The vector equation (1) can be considered as a system of n linear 
equations with n + 1 unknowns; some or all of the entries of the 
noise vector, and the gain. The size n of the system depends only 
on the number of non- vanishing coefficients in the original signal 
and the chosen basis representation (through W(l)), and the co- 
efficients on G is directly controllable via s 0 . This means that the 
linear equation system can be tailored to fit an approximate solu- 
tion method such as least square. Tf the noise is normal distributed, 
e t ~~ Nfyi, a 2 ), this will give the best result, independently of the 
wavelet transform. Inverse transforming (1) yields 



s r = Gst + e t . 



(2) 



Note that s t can be any signal, and s r is the result of the trans- 
mission. A least square approach could be formulated through a 
rewriting of (2) to 

\\*r-Gs x - »l\\ 2 =a*N t (3) 

where N is the length of the signal. Let a be the smallest value for 
which (3) holds. Then (3) is an elliptic paraboloid with minimum 
in the (G, fi) plane. This minimum point is found when the G and 
M discriminants are zero simultaneously. Solving that yields 



G = 



(s r ,s t ) - /x<s t ,l) 



with 



(4) 



(5) 



(6) 



(s r ,s t )(s t ,l)-f ||s t |l 2 ( Sr ,l) 
(s t ,l> 2 -7V-||s t |p 

The smallest a is then 

, tf ( St ,l) 2 -(s r ,s t ) a + |ls t || 3 Hs f H 2 a 

IMP* 

As long as the noise is normally distributed these estimations are 
fairly accurate. Unfortunately there is no way of immediately 
telling how accurate the estimation of G is. 

Another way is to exploit two properties of the linear equa- 
tions, namely that G is in all them and that the noise, although 
unknown sample by sample, has some approximately known prop- 
erties (mean, energy, frequency distribution etc.). With the wavelet 
transform it is then possible to make the G estimation more robust 
than in the previous solution; having the freedom in choice of orig- 
inal signal, another possibility is to design the original signal such 
that (e w t , so) is close to zero. Then 



(s wr >so) _ (e w t, so) ^ (s Wf ,sp) 

(S„,S 0 ) <S„,S 0 ) ~~ {s n ,S0> ' 



(7) 



\ 



) 



If certain frequency bands in the noise are 'nice', i.e. has low mean 
square error or contains nearly white noise, the wavelet transform 
can target these intervals making, it easy to design a signal s 0 
which has high probability of (e wt ,s 0 > being close to zero. 

2.3. Verification of measurements 

The special structure of the systems allows for solutions involv- 
ing only easily calculated quantities, like the inner products. Such 
a simple approach, however, is sensitive to time localized noise, 
spikes in particular. To reveal occurrences of this type of noise the 
freedom in design of the original signal comes in handy. Let s 0 
through sk be K -f 1 signals with the properties 

(s*,s m )=0 fork^m and <s fc ,ewt) % 0, 

where e t is a typical noise occurrences. When transmitting s m the 
other K signals, called verifiers, can be used to indicate the quality 
of an estimation on the form (7). A time localized noise occurrence 
is still time localized in a wavelet transform, and <s wr , s m ) could 
in such a case very well be a poor estimate of the gain. The inner 
products (s wn sjk> for k ^ m would be similarly affected by the 
noise, independently of the channel gain. Comparing <s wr ,s m ) 
with each of the verifiers inner product with s wr gives a quality of 
the estimation of G. A comparison could be done like 



p = 1 - exp 



/ -BK(s WT ,s m ) 2 ' 

K 

( S wr,Sjfc) 2 



(8) 



where the exp is introduced to make the quality of the estimation 
absolute, that is 0 < p < 1. Whenever the verifiers detect en- 
ergy in the received signal (energy which then have to originate in 
noise), p becomes small. The higher the ratio of estimated gain and 
verifier energy, the closer p gets to 1, with B controlling how close. 
When some parts of the received signal is corrupted by noise, there 
might still be other usable parts. Then (7) and (8) can be restricted 
to these parts to give a G estimate unaffected by the corrupting 



3. RESULTS 

A number of suggestions on how to utilize the wavelet transform 
to measure channel gain has been given in the previous section. 
The result of using a few of them in an infra red emitter-receiver 
implementation is presented in this section. In all cases the result 
of wavelet transforming with W is the fourth level in a wavelet 
packet transform decomposition. This level consists of eight ele- 
ments, as indicated in the figures. The filler is Symlet 8 [2], and 
due to the short filter periodization is used. The noise is recorded 
with the infra red photodiode at 5 kHz, and the noise vector e t is 
shown in figure lc. As original signal a step function is chosen, 
figure la, and the function x lOOOx + 2048 is chosen as affine 
mapping. This mapping applied to s 0 is shown in figure lb, which 
is then the emitted signal. The gain is known to be 0.00642. The 
received signal in shown in figure Id, and the transform of this is 
shown in le. 

Applying the least square approach (4)-(6) to the transmitted 
signal s t , see figure lb, and received signal, see figure Id. gives 
a gain and mean which deviates only very little from the true val- 
ues (see table I). However, the addition of a single spike in the 




0 64 128 192 256 320 384 448 512 



1360 



0 64 128 192 256 320 384 448 512 




128 192 256 320 384 448 512 

Figure 1 : Experimental data from an infra red emitter-receiver im- 
plementation, a) original signal, b) emitted signal, c) transmission 
noise (experimental data), d) received signal, e) the received signal 
wavelet transformed. 



received signal (in this case 200 is added to sample 177) makes 
the estimation of G rather poor, and there is no immediate way of 
detecting this. 

The detection of a poor estimate is possible with the inner 
product approach (7). First, applying it to the wavelet transform 
of the received signal (figure le) yields a good estimation of G 
(see table 2). To determine the quality of the estimation two ver- 
ifiers, sj and s 2 in figure 2a and b, are used. They have the same 
structure as s 0 while at the same time being orthogonal to s 0 . The 
absolute quality measure p in (8) also shows a good estimation (£? 
is chosen such that p = 1 - exp(-B ■ 200) = 0.9). When applied 
to the transform of the spiked, received signal (the same spike as 
before), the estimation is poor, which is easy to tell from the p 



Description 


Fig 








No extra disturbance 


Id 


0.2% 


0.03% 


5% 


Added 200 to sample 177 


2d 


23.3% 


0.29% 


171% 



Table 1: Results of applying the least square methods (4)-{6). The 
three percentage columns shows the deviation between the esti- 
mated and the real values. 




Table 2: Results of applying the solulion method (7) The ihirH 
£*,.„). and the last column the corresponding quality measure 
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If a number of spikes occur distributed along the entire re 
caved s.gn*. the above method might fail. In this case a sp ke « 
duction method ,s useful. The estimated value of G (^ich fw ,hU 
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4. CONCLUSION 

oases, implementations etc., and the freedom in r . . 
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